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Abstract 



An improved quasistatic model is used to describe the ionization process of 
atoms in intense linearly polarized fields. Numerical calculations of Above- 
Threshold-Ionization (ATI) energy spectra and photoelectron angular distri- 
butions (PAD) of hydrogen atoms reveal clearly plateau and side lobe struc- 
tures, which are in good agreement with recent ATI experiments. Our results 
show that, the existence of this plateau and side lobes is consequences of the 
classical kinematics of the electrons in combined atom and laser helds. Fur- 
thermore, we find that the onset of these unusual phenomena is related to the 
onset of the chaotic behavior in the system. 
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I. INTRODUCTION 



The most recent important finding in Above-Threshold-Ionization (ATI) experiments is 
the existence of a plateau formed by the high-order ATI peaks which halts and even tem- 
porarily reverses the decrease of the heights of the peaks with increasing order [1]. Moreover, 
the Photoelectron Angular Distribution (PAD) in the transition region at the onset of the 
plateau, seem unusual: while the angular distributions both below and well above this region 
are strongly concentrated in the field direction, additional 'side lobes' at angles between n/6 
and 7r/4 with respect to the field direction develop in the transition region [1]. 

It has been concluded by numerical calculation of the one-electron Schroedinger equation 
that the new phenomena can be understood in the context of single- electron ionization 
dynamics [1]. The conjecture that the electron returning to the ion core will cause a number 
of observable consequences , such as the cut-off law of the high-order harmonic production, 
has been made by many authors [2, 3] and it now seems evident that recent findings of 
ATI result from the rescattering effect as well. A fully quantum mechanical calculations 
for the rescattering effect based on a two-step model with a simplified delta source has 
been well performed by Becker et al [4]. An extensive analysis of this effect based on the 
Keldysh-Faisal-Reiss (KFR) theory with realistic atomic potential has also been presented 
[5]. Lewenstein and coworkers have made analysis on lobes in ATI and its relation to 
rescattering by using semiclassical method [6]. However, the precise physical origin of the 
plateau and lobes is still unclear. Are they due to some genuinely quantum mechanical 
effect such as interference or can they be traced back to the essentially classical behavior of 
the electrons in the combined atom and laser fields? This question is partially answered by 
Paulus et al [7]. However the emission of photoelectrons in the direction of the laser electric 
fields is strongly underestimated by their simple classical model. 

In this paper, first we shall develope a 3D quasistatic model which generalizes the well- 
known quasistatic model [2, 3] by including the effect of Coulomb potential on electron motion 
after tunnel-ionization. Taking the simple Hydrogen atom as an example, we shall investigate 
its ionization process in an intense laser field based on the improved model. Our results 
show clearly that, the plateau structure in ATI spectra and side lobes in the PADs can be 
attributed to the pure classical behavior of electrons in the combined Coulomb and laser 
field. Second, we shall present numerical evidence that, the unusual phenomena observed 
in ATI data result from a kind of irregular motion in the nonlinear dynamical system. Our 
discussions connect the optical field ionization problem with another active field, namely, 
Chaotic Dynamics. This connection will enrich and deepen our understandings on atom-laser 
interaction. 

II. AN IMPROVED QUASISTATIC MODEL 

As Keldysh parameter 7 = (Io/21/p) 1 ^ 2 <C l,Tiu/2U p <C 1 and field strength F < F th , 
tunneling ionization occurs. Here, h is the Planck constant, Iq is the ionization potential of 
the electron, uj is the field frequency, Z the core charge, U p = \^ ^ 2 the ponderomotive po- 
tential, F t h the threshold field (its meaning and value will be given in following disscussions) . 
The probability of the tunnel-ionization is first obtained by Landau [8] for a hydrogen atom ( 
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in the ground state) and then be extended to complex atoms [9]. More recently, Delone and 
Krainov derived an analytical expression for the probabilities of tunnel-ionization of atoms 
and for the energy and angular electron spectra in a strong low-frequency electric field [10]. 
Their discussions give a perfect description of the tunneling problem and will be helpful to 
provide the initial conditions of an electron after tunneling in our 3D model. 

In the second step of the quasistatic procedure, by including the Coulomb potential 
as well as the laser fields, we obtain a complete Newtonian equation, which describes the 
motion of an electron after tunneling ionization. 

Ze 2 f 

+ (—ZeFcosujtel + (3ZeFsinutay) , (1) 



where /3 — for linearly polarized laser and (3 — ±1 for circular polarization. 

It is convenient to introduce the compensated energy E c advocated by Leopold and 
Percivalfll] (the E c is just the Hamiltonian function in the velocity gauge), 

E c = — -[f + (ZeF/m e u) sm(u>t)e z - (3(ZeF/m e u) cos(iut)e y ] 2 - Ze 2 /r . (2) 

When an electron is ionized completely, the Coulomb potential is weak enough and E c tends 
to be a positive constant which is just the ATI energy (Eati) m ultrashort pulse laser. 

In what follows, taking the simple Hydrogen atom in the linearly polarized fields as an 
example we shall demonstrate this approach. For simplicity, the atomic units are used ( that 
is, m e — e — h — la.u ). The field parameters are chosen as F = 0. 06a. u. (1.2 x 10 14 w/cm 2 ) 
and tu = 0.04242(A = 1064nm), I = U p = 0.5a.«.(13.5eV), then 7 = 0.71 < 1 and 
hu/2Up = 0.04242 << 1, and therefore this situation is well in the tunnel ionization regime. 

We first recall the deduction of the effective potential in [8] for the self-consistence of the 
paper. Schroedinger's equation for the hydrogen atom in a uniform electric field is of the 
form 

dA + I + --ez)^ = 0, (3) 

where A is the Laplace operator; e = Fcos^o, the quasistatic field; 0o = wto, the field phase 
at time to- 

Let us introduce the parabolic coordinates £, 77, ip 



x = V^cos^,!/ = y^sin^z = (f - rj)/2. (4) 
We seek the eigenfunctions in the form 



ft 

where m is the magnetic quantum number. 
Then equation (3) is brought into the form 



d 2 fi ( h di m 2 - 1 1 

dp" + ( 2- + T-^ 4^/1 = 0, (6) 
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d 2 f 2 , Jo , fo m 2 -l 1 

/3i + & = l. (8) 

We shall regard the energy I as a parameter which has a definitive value. Solving the 
equations, as functions of Jo and e, and then the condition (3± + 02 — 1 will give the energy 
as a function of the external field. 

Above equations show that there exist a potential barrier along the rj coordinate, the 
ionization of the electron from the atom in the direction z — > — oo corresponds to its pas- 
sage into the region of large rj. To determine the ionization probability it is necessary to 
investigate the form of the wave function for large rj and small £. Neglecting the stark shift 
we have I = — \ and (5\ — 02 — \ and m — 0. In the absence of the field, the wave function 
is 

^ = e-K+W 2 /^. (9) 

When the field is present, the dependence of the \I/ on £ can be regarded as being the same 
as (9), while to determine its dependence on r\ we have the equation (from (7)), 

d 2 ® , 1 1 1 1 ^ 

77T + ~7 + 7~ + 7^ + 7 er ? $ = ' 10 

where $ = y^^. 

Above equation has the form of the one-dimensional Schrodinger's equation with the 
potential U(r)) = — ^ — — |er] and the energy K = — |. The turning point ( at the 
outer edge of the suppressed potential ) where an electron born at time to is determined by 
U(rj) = K and expressed by 



l - + 2s-^ cos 6, s = *njfr, 6 = \ arccos(-f ), 



770 3e ' \ ' :5 ' v 2s' 

1 2 2 2 1 

As e > F t ^, the turning point will be complex, which determines the threshold value of the 
field F th = 0.113a.u.. 

WKB method is used to obtain the tunneling rate of the above system and gives [8] 
w(0) = ^exp(—-^), which is used in the original quasistatic model[2,3] (quasi one- 
dimensional model). This model can qualitatively predict the ATI spectrum and provide a 
analytic expression of the initial parameters in the problem of the optical field ionization 
X-ray laser. However, this simple model can not include the rescattering effects that result 
in many important observable phenomena, such as Plateau in ATI spectra. Morever, it can 
not predict the PADs for it is an one-dimensional model. Recently, Delone et al generalized 
the above tunneling formula and obtained the probabilities of tunnel ionization for both 
energy and angular energy as follows[10], 
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w(p±) = w(0)w(l),w(l) = exp(—p 2 ± /e) 



(12) 



where p± is the perpendicular momentum of the photoelectrons. 

The above formula will be used to weight an classical trajectory in our three-dimensional 
quasistatic model. 

Since the system is azimuthal symmetric about the polarization axis, we can restrict 
the motion of an electron on the plane (x, z). In addition, for the exponential decay of the 
tunneling rate along the £ axis, it is a reasonable assumption that £ — for an electron 
just after tunneling. Thus, the initial position of an electron born at time t is given by 
xo = 0, z = —r]o/2 from equation (4). The initial velocity is set to be v z — 0, v x — v x o- The 
weight of each trajectory is evaluated by the Delone's expression of tunnel rate (12), 

M0o,^o) = w(0)w(l),w(0) = -exp(~) f w(l) = —exp(-v 2 x0 /e) . (13) 

e oe en 

In applying Delone's formula to our problem, the second term should be normalized to 
present form. From this formula, one can readily verify that, for a fixed initial perpendular 
velocity v x q, the trajectory corresponding to the initial field phase <po — (i.e. at maximum 
quasistatic field e = F) has the maximum probability; for a fixed initial field phase <po, the 
trajectories corresponding to v x q = ^/e/2 has the maximum probability. 

In our computations, 10 5 initial points are randomly distributed in the parameter plane 
— 7r/2 < o < tt/2,v x0 > so that the weight of the chosen trajectory is larger than 1CT 11 . 
Each trajectory is traced for such a long time that the electron is actually ionized (this can 
be evaluated by using the compensated energy Eati-) The final ATI spectra and the PADs 
can be obtained by making statistics on an ensemble of classical trajectories. The results 
have been tested for numerical convergence by increasing the number of trajectories. 

Before we turn to our numerical results, we would like to mention that similar model has 
been used successifully to double ionization of He atoms by Brabec et al [12]. Differently, 
in this model, the initial conditions of electrons after tunneling is derived strictly from 
the effective potential in Landau's tunneling theory rather than Coulomb potential. This 
deduction gives the correct threshold field F t h = 0.113a.u for 3d hydrogen atom. This model 
is restricted to the situation where the external field strength is smaller than the threshold 
value. 



III. RESULTS 

A. ATI Spectra and PADs of Photoelectrons 

In Fig.l, we demonstrate the ATI spectra and the total angular distribution (of the 
emitted electrons with respect to the angle 9, i.e. at the detector) calculated from our 
model. The original simple quasistatic model predicts an sharply decreasing ATI spectra 
curve and the photoelectron energy can not be larger than 2C/ p [2, 3]. Our results show 
clearly that the rescattering increases the fraction of the hot electrons. This is due to that 
an electron has a higher probability of staying in the vicinity of the nucleus and then absorb 
more photons in the rescattering process. In particular, the ATI spectrum exhibits a sharply 
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decreasing slope (region I, - 2U P ) followed by a plateau (the transition region II, 2U P - 8U P ) 
and again a sharply decreasing slope (region III, > 8U P ). The height of the plateau is three 
orders of magnitude below the maximum of the spectrum and its width is about 6U P . This 
phenomenon is qualitatively in agreement with the existing experiments and various theories 
[1, 4, 5, 6]. The distinctive feature is the plateau's extension about QU P width to rather high 
energy before abruptly delcining at about 8U P . This differs from the previous experiments 
[1] but is much closer to recent experiment [13] which is performed in strong tunneling limit. 
The total angular distribution in Fig. lb contains only the electron initiated in phase interval 
[— 7r/2, 7r/2]. The angular distribution of the electrons originated in [7r/2,37r/2] will be the 
mirror image with respect to 90°, so the sum of the two contributions will show a main 
concentration in the field direction. 

Furthermore, in Fig.2 we calculate statistics on the angular distribution of photoelectrons 
in three different energy regions. The most striking feature of the plots is the existence of 
a slight slope up to 40° followed by a sharp cut-off, i.e. no photoelectrons in the transi- 
tion region emit at angles much larger than 40° (see fig. 2b). This remarkable phonomenon 
corroborates the data of Paulus et al [14] , which we think is due to the pure tunneling 
nature. In contrast, in the mixing regime where multiphoton ionization becomes significant, 
the angular distributions show no such cut-off, and there appears to be emission even at 
90°. Considering the rescattering effects from a simple classical model, Paulus et al find a 
peak at 30° of the angular distribution for the electrons in the transition region. However, 
the emission of the photoelectrons in the direction of the laser electric field is largely under- 
estimated [7]. In our model, a complete Newtonian equation is used to simulate the electron 
motion after the tunneling so that the effects of multiple returns are included. Therefore, 
our model gives more reasonable results on the emission of the photoelectron in field di- 
rection. Furthermore, the PADs in the regime of plateau exihibit additional peak which 
is so called side lobe phenomenon observed in the experiments. One example is shown in 
figure 2d for 4U P ± 0.1U P , there the side lobe occurs at about 40°. The sum of the side 
lobes for photoelectron in the plateau region are responsible for the slight slope followed by 
a sharp cut-off observed in the Fig. 2b. Our numerical results also show that the angular 
distributions both below and well above the plateau region are strongly concentrated in the 
field direction. Detailed investigations show that, aside from these electrons which directly 
drift away without returning to the core , the most electrons in the region I experience 
the forward scattering , while the electrons in the region III are backscattered by almost 
180°. However, in the plateau region the forward scattering and back scattering are equally 
important and seem have equal probability. 

In Fig. 3 we demonstrate three typical trajectories of the photoelectrons corresponding 
to the region I, II and III in Fig. la respectively. Generally speaking, after tunneling, most 
electrons will be driven by the external fields to return to and interact with the core. This 
rescattering process greatly determine the final energy and momentum distributions of pho- 
toelectrons. Without considering the rescattering, the PADs are expected to become more 
peaked along the polarization axis as the order of the ATI peak increase, and the ATI 
spectrum be a sharply decreasing curve. However, in the rescattering process, an electron's 
energy changes a lot through an exchange of momentum with the core. It can obtain high en- 
ergy from a strong collision. Then, a near 180° backscattering which imply a large exchange 
of momentum certainly gives the highest energy (region III). A forward scattering with small 
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emission angle only provides a relative low energy for an electron (region I). Whereas, in 
the transition region ( region II) the situation become quite complicated. Multiple return 
and long-time trapping can be experienced by those electrons in this region. The classical 
trajectories show complex behavior. This fact leads to the randomness in emission direction 
of photoelectrons. As will be seen later, this is a signature of the chaotic behavior of the 
trajectories that constitute the plateau region.. 

B. Chaotic Behavior in Rescattering Process 

As was shown above, in our model, the rescattering process of the electrons after tun- 
neling can be well described by the Newtonian equation (1). ATI spectra and PADs can be 
obtained by calculating statistics on an ensemble of trajectories correspinding to different 
initial field phase and perpendular velocity. To investigate the detailed dynamical mecha- 
nism underlying those unusual phenomena such as plateau structure in ATI data, we shall 
first fix the velocity and calculate the initial phase (0o) dependence ATI energy and emission 
angle (Fig.4,5). As the initial velocity perpendicular to the polarization of the electric field is 
relatively large (case a and b), we have smooth phase dependence ATI energy and emission 
angle which imply regular motion of classical trajectories. In these cases, an ionized electron 
with lower ATI energy tends to possess a higher probability (near O — 0), and the plots 
of angular distribution indicate that hoter photoelectron tends to emit in the polarization 
direction. Moreover, in these cases the ATI energy of the photoelectron is smaller than 2U P 
approximately. Considering the exponential form of the expression of ionization rate (13), 
one can conclude that these cases correspond to a sharply decreasing spectrum such as the 
region I in Fig. la. Therefore these trajectories have very small contribution to the unusual 
distributions in the plateau region in Fig. la. For case c and b , things become quite different. 
The dependence of Eati and the emission angle on the initial phase is poorly resolved in 
a region near zero point where many strong peaks are observed. Successive magnifications 
(Fig. 6) of the unresolved regions show that any arbitrarily small change in the initial phase 
may result in a substantial change of the final electron energy and emission angle. Multiple 
returns and even infinite long time trapping can occur in this region. This is the character 
of chaotic behavior, namely, sensitive dependence on the initial condition, which has been 
observed in many chaotic scattering models [15]. The usual rule that the higher energy cor- 
responds to the lower emission angle and the lower probability which is identified in the case 
a and b is broken up now. In the unresolved regions, some chaotic trajectories can produce 
rather high ATI energy (> 2U P ) and therefore are responsible for the unusual structure - 
plateau in ATI spectra. Additionally, emission angle of these trajectories are also strange, 
this fact leads to the unusual angular distribution in the transition region. In the unre- 
solved region, the electrons with chaotic trajectories are backscattered or forward scattered 
randomly, this also explains the fact observed in the last section , that the coexistence of 
forward scattering and back scattering motion of the photoelectrons in the plateau region. 

Now we turn to apply the nonlinear dynamics theory[16] to give a quantitative description 
of the singularity of the unresolved regions in the phase dependence ATI energy plot. We 
call a value 0o a singularity if, in any small neighborhood domain there is a pair of <j>o 
which have different signs in their final moments projected in the field direction, i.e. they 
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correspond to the forward scattering or back scattering, respectively. This singularity results 
from the unstable manifolds and stable manifolds which intersect with each other and form 
a nonattracting hyperbolic invariant set with fractal structure. Taking the situation of 
Fig. 4c as an example, we employ the uncertainty exponent technique to obtain the fractal 
dimension of the singular set in the phase axis. We randomly choose many values of 0o 
in an interval containing the fractal set. We then perturb each value by an amount e and 
determine whether the final canonical momentum projected in the polarization direction 
corresponding to initial (fi ,(fio — e and <fio + e have the same sign. If so , we say that the 
(f) value is e-certain; if not, we say it is e-uncertain. We do this for several e and plot on 
a log- log scale the fraction of uncertain O values /(e). The result is plotted in fig. 7 which 
shows a good straight line and indicates a power law dependence /(e) ~ e 7 , where 7 = 0.25. 
The exponent 7 is related to the dimension of the fractal set of the singular <p s values in 
0o-axis in Fig. 4c by [16] 

D = 1 - 7 = 0.75 . (14) 

As discussed above, the irregular phase dependence ATI energy and emission angle result- 
ing in the unusual phenomena such as plateau and side lobe in recorded ATI data, come from 
the chaotic scattering behavior in system (1). To demonstrate phase structure of the chaotic 
scattering, we choose the initial condition as those points in phase space well in the asymp- 
totic region of the potential, typically, x = 0,v x G (—0.5,0.5) and z = 150, v z G (—2,0) and 
then trace the evolution of the 5000 trajectories. The stroboscopic map for such an ensemble 
of scattering trajectories is shown in Fig. 8. Since the system can flow arbitrarily between 
the bounded and unbounded regimes, it is a non-compact system. The phase sections show 
a rich and complicated dynamical structures which confirms the existence of the chaos. 

IV. CONCLUSIONS AND DISCUSSIONS 

In this paper, we calculate the ATI spectra and angular distribution in the regime of 
tunnel-ionization from an improved quasistatic model , and show that those unusual phe- 
nomena recorded in ATI data can be traced back to the essential classical behavior of the 
electrons in the combined atom and laser fields. The properties of ATI spectra and angu- 
lar distribution presented in this paper reflect some fundamental features of ionization in 
the strong field tunneling limit. For example, recent high precision measurements of helium 
photoelectron energy and angular distributions in the pure tunnel regime demonstrate much 
similar properties, such as the flat energy distribution that extends out to surprisingly high 
energies before abruptly truncating at 8 — 10C/ p [13]. In another aspect we associate the un- 
usual phenomena in ATI spectra and PADs in a real atom system with the irregular chaotic 
behavior. We think that nonlinear dynamics theory is very helpful in understanding the 
behavior of atoms in the intense laser fields just like its successful applications in photoab- 
sorption spectrum of Rydberg atoms [17] and the microwave ionization of the highly excited 
atoms [18] 
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CAPTIONS OF FIGURES 

• Fig.l ATI spectra and total PADs calculated from our model. F = 0.0Qa.u.,uj = 
0.04242a.w.. 

• Fig. 2 PADs for different energy regions. a,b and c correspond to the region I, II and 
III in fig. la respectively. The d shows the PAD at energy around AU P , where the side 
lobe near 40° is very obvious. 

• Fig. 3 Three typical trajectaries correspond to the region I, II and III in fig. la re- 
spectively. a)Initial conditions are t = 8. 050a. w., v x0 = 0. 123a.it., the weight 
2.839 x 10~ 4 a.-u. The final ATI energy is 0.29a.w., emission angle —3.2°. A weak 
collision occurs in this case. b)Initial conditions are t = —6.356a.u.,v x0 = 0.130a. it., 
the weight 3.645 x 10" 4 a.w.. The final ATI energy is 3. 58a. it., emission angle —13.3°. 
Multiple return occurs in this case. c)Initial conditions are to = 0.871a.u.,v x o = 
6.00 x 10~ 5 a.w., the weight 3.154 x 10~ 7 a.w.. The final ATI energy is 5.02a.w., emis- 
sion angle —178°. A strong collision occurs in this case. 

• Fig. 4 Phase dependence ATI energy spectra for four different initial perpendicular 
velocities. F = 0.06a.w., cu = 0.04242a.w.. 
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Fig. 5 Phase dependence PADs for four different initial perpendicular velocity. F = 
0.06a.M.,o; = 0.04242a.w.. 

Fig. 6 Successive magnifications of fig4c. 

Fig. 7 Scaling Law of the Singular Set. 

Fig. 8 Stroboscopic map a) in plane (x,v x ), b) in plane (z,v z ) for 5000 trajectories 
originating asymptotically at x = 0,v x G (—0.5,0.5) and z = 150, v z G (—2,0). 
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FIG. 1. ATI spectra and total PADs calculated from our model. F = 0.0Qa.u.,uj = 0.04242a.u.. 
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FIG. 2. PADs for different energy regions. a,b and c correspond to the region I, II and III 
fig. la respectively. The d shows the PAD at energy around 4U P , where the side lobe near 40° 
very obvious. 
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FIG. 3. Three typical trajectaries correspond to the region I, II and III in fig. la respectively. 
a)Initial conditions are to = 8.050a.u.,v X Q = 0.123a.u., the weight 2.839 x 10~ 4 a.u. The final 
ATI energy is 0.29a.u., emission angle —3.2°. A weak collision occurs in this case. b)Initial 
conditions are to = — 6. 356a. u., v x q = 0.130a.u., the weight 3.645 x 10 -4 a.u.. The final ATI energy 
is 3.58a.u., emission angle —13.3°. Multiple return occurs in this case. c)Initial conditions are 
to = 0.871a.u., v x0 = 6.00 x 10' 5 a.u., the weight 3.154 x W 7 a.u.. The final ATI energy is 5.02a.u., 
emission angle —178°. A strong collision occurs in this case. 
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FIG. 4. Phase dependence ATI energy spectra for four different initial perpendicular velocities. 
= 0Ma.u.,LO = 0.04242a.u.. 
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FIG. 5. Phase dependence PADs for four different initial perpendicular velocities. 
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FIG. 7. Scaling Law of the Singular Set. 
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FIG. 8. Stroboscopic map ci) in plane (x, v x ), b) in plane (z, v z ) for 5000 trajectories originating 
asymptotically at X = 0,v x € (-0.5,0.5) and Z = 150, v z £ (-2,0). 
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